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• ABSTRACT 

00 ■ 

^-H Aims. We propose a new approximation for the cyclo-synchrotron emissivity of a single electron. In the second part of this work, we discuss a 
simple application for our approximation, and investigate the heating of electrons through the self-absorption process. Finally, we investigate 
the self-absorbed part of the spectrum produced by a power-law population of electrons. 

Methods. In comparison to earlier approximations, our formula provides a few significant advantages. Integration of the emissivity over the 
whole frequency range, starting from the proper minimal emitting frequency, gives the correct cooling rate for any energy particle. Further, 
the spectrum of the emission is well approximated over the whole frequency range, even for relatively low particle energies (/3 O.i), where 
O I most of the power is emitted in the first harmonic. In order to test our continuous approximation, we compare it with a recently derived 
' approximation of the first ten harmonics. Finally, our formula connects relatively smooth to the synchrotron emission at /? = 0.9. 
^ I Results. We show that the self-absorption is a very efficient heating mechanism for low energy particles, independent of the shape of the 
■ particle distribution responsible for the self-absorbed synchrotron emission. We find that the energy gains for low energy particles are always 
, higher than energy losses by cyclo-synchrotron emission. We show also that the spectral index of the self-absorbed part of the spectrum at 
• • very low frequencies differs significantly from the well known standard relation I(v) ~ v''^. 
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1. Introduction 

Synchrotron emission is well understood (see e.g. reviews by 
Ginzburg & Syrovatskii 1965, 1969; Pacholczyk 1970) and is 
thought to be responsible for a significant part of the radiation 
we receive from a variety of cosmic objects, such as super- 
nova remnants, radio jets, compact radio sources, active galac- 
tic nuclei, and gamma-ray bursts. The reverse process, syn- 
chrotron absorption, has only recently disclosed some novel 
features, when the attention of researchers shifted from what 
is the amount of absorbed flux of photons to what happens 
to the absorbing electrons. In fact, they can absorb the en- 
ergy of the photons and thereby change their initial distribu- 
tion (GhiseUini, Guilbert & Svensson 1988, hereafter GGS88), 
and/or absorb the momentum of the photons, with the possibil- 
ity of bulk motion acceleration (Ghisellini et al. .l990j . 

These authors demonstrate that the only stable equilibrium 
solution of particles emitting and absorbing synchrotron radi- 
ation is a relativistic or quasi-relativistic Maxwellian distri- 
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bution. This paper ended a long debate about the existence 
of another equilibrium solution: a power law of slope 3, i.e. 
N{y) oc y^^, which was the main result of the so-called "Plasma 
Turbulent Reactor" (PTR) theory, as introduced in a series of 
papers in the 1970s (Norman I 19771 Norman and ter Haar l 19751 
Kaplan & Tsvtovich ll973> . Note also that the stability of this 
N{y) oc y^^ solution was already questioned by Rees (1967), 
stating that this power law solution would evolve away from 
y^^, if slightly perturbed (see also the numerical results by 
McCray 1969, demonstrating this instability). 

One of the aims of this paper is to explicitly demonstrate, 
that the y^^ distribution is not only unstable, but is not even an 
equilibrium solution. To do so in an accurate way, however, 
it is necessary to also consider the trans-relativistic regime, 
namely, the cyclo-synchrotron emissivity and absorption co- 
efficient. This, of course, is the more complex regime, because 
the emitted power is not concentrated at all in the first harmonic 
at the typical synchrotron frequency (i.e. ~ y^VB, where vb is 
the Larmor frequency). Recently, a useful approximation has 
been proposed by Marcowith & Malzac (2003), introducing 
polynomial expressions for the first 10 harmonics for a range 
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of particle energies. They compare their results with an exist- 
ing analytical formula that tried to approximate the emission 
(and the absorption) with smooth functions (i.e. not as sums 
of harmonics), as proposed by Ghisellini, Haardt, & Svensson 
(1998). From this comparison it appears that there might still be 
room for improvement in this smooth, approximated function, 
which is the second main aim of our paper here. 

We present our new approximation in Sect. 2, and com- 
pare it with the Marcowith & Malzac (2003) results. We show 
that our approximation works well for particles with /3 < 0.9, 
where f5c is the particle velocity, corresponding to y < 2.3. 
For slightly higher energies, the standard synchrotron formu- 
lae describe the shape of the emission well, but the frequency- 
integrated emissivity must still be corrected to become equal to 
the cooling rate (which is known exactly for any particle en- 
ergy). Since we are interested in the total amount of energy 
absorbed and lost by a single particle, in Sect. 3 we introduce a 
correction to the standard synchrotron formula, which is impor- 
tant for 'y>2.3, but automatically negligible for ultra-relativistic 
energies. 

We then consider the rate of energy gains and losses suf- 
fered by an electron of a given energy as a consequence of 
cyclo-synchrotron emission and absorption (Sect. 4), show- 
ing that only particles at a single energy can be in equilibrium, 
where gains equal losses, independently of the slope of the par- 
ticle distribution that produces the cyclo-synchrotron intensity. 
Studying the low frequency part of the synchrotron intensity in 
detail, we show that there are novel features below the Larmor 
frequencies that appear to have been overlooked in the past. 
These are presented in Sect. 5. Finally, we draw our conclu- 
sions in Sect. 6. 



2. Approximation of the cyclo-synchrotron power 
spectrum 

The single particle cyclo-synchrotron power spectrum can be 
approximated relatively well by a simple analytical formula. 
One of the best approximations was proposed by Ghisellini, 
Haardt, & Svensson ( fT998t 



Pc{y,p) 
fip) 



4 0-tCUb 2r, N 

3 Vb 



fip){l-- 

VB 



1 + ap^ 



(1) 
(2) 



where a = 3, vg - eBKlnm^c) is the Larmor frequency, 
Ub - B^/i^n) the magnetic field energy density, B the mag- 
netic field strength, p - y/3 - yjy^ - 1 the dimensionless par- 
ticle momentum, y the particle Lorentz factor related to the to- 
tal energy hy E - ynieC^, /3 the particle velocity in units of c, 
and CTj, ifie, c are constant Thomson cross-section, electron 
rest mass, and the velocity of light, respectively. Eq 1 describes 
the power, integrated over the emission angles, of an "average" 
electron: i.e. the power has been averaged over the pitch angle, 
which is assumed to be distributed isotropically. In this case 
the emissivity of a single electron is equal to the emitted power 
divided by the solid angle factor An. We use the term "power 
spectrum" to indicate the emitted power as a function of fre- 



quency of an "average" electron, in the sense specified above. 
This phenomenological formula has three advantages: 

+ can be integrated easily over the frequency range 
-I- integration from v = up to infinity gives the correct cool- 
ing rate 



ycirieC^ - -o-jcUbP^, 



(3) 



+ has the correct frequency dependence [exp(-v/vB)] at large 
harmonics (v » vb) 

On the other hand, the formula also has three significant disad- 
vantages: 

- gives the correct cooling rate only if integrated from v - 
Vb, therefore cannot correctly describe the emission level 
below the frequency vb 

- does not approximate the emission spectrum well for /3 < 
0.5, as we will show 

- does not join smoothly to the synchrotron power spectrum. 

Therefore, in order to improve this formula, we introduce three 
important modifications. 

First, the problem with the lower integration limit, (v - vb) 
is solved by multiplying the formula with the term 



g(y,p) = 



V - Vmin(p) 



(4) 



where Vmin indicates the minimal emitting frequency of the first 
harmonic 
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This modification significantly reduces the power emitted be- 
low the frequency vb in comparison to the original formula and 
provides an automatic cut-off at the limiting frequency (vmin). 
Note that this additional term becomes unity for v » v^in, thus 
maintaining the original formula in this regime. However, we 
introduce a v-dependent term in front of the exponential func- 
tion (exp[/(l - v/vb)]), which also depends on the frequency. 
Therefore, the integral over the frequency range of the new for- 
mula must be expressed by the exponential integral. 

Second, to improve the shape of the spectrum for /3 < 0.5, 
we replace the term f{p) by a modified expression f'(p) read- 



mg as 
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where a = 3.65 and b = 0.02. The term (p^ + b)/p^, which 
makes the diff'erence between f{p) and f'ip), becomes equal 
to unity for p '» b, and therefore our approximation becomes 
equivalent to the original formula in this regime. The necessary 
changes to the constant a = 3 — > 3.65 is discussed in the next 
section. 

Third, the new expression is normalized in order to yield 
the correct cooling rate (Eq. |5} when integrating over the fre- 
quency range. This normalization is done by multiplying with 
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a factor 
c{p) = ■iexp 



vb 
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where Ei is the exponential integral (e.g. Press et al. ll989> . 

Finally, the improved approximation for the cyclo- 
synchrotron power spectrum of a single particle is given by 



Pciv,p)^- 



P c(p)g{v,p)f (p)exp 



3 Vb 



(8) 



In comparison to the original relation this formula has three 
significant advantages: 

+ integrated from v - v^in gives the correct cooling rate for 

any energy particle 
+ describes the emission spectrum well, starting from the 

minimal frequency for a wide energy range (J3 - 0.01 — > 

0.9) 

+ provides a relatively smooth connection to the synchrotron 
power spectrum at /? = 0.9. 

The price we pay for these improvements is that the inte- 
gral over frequency of the new formula must be expressed in 
terms of an exponential integral. Moreover, in some sense our 
formula averages over the harmonics, providing a continuous 
emission spectrum. Therefore, it cannot be used for model- 
ing the emission where the cyclo-synchrotron lines are ob- 
served directly (e.g. Pottschmidt et al. 2005'). On the other 
hand, our approximation may have a wide range of applica- 
tions in sources where the particle energy extends at least over 
one order of magnitude. In such a case the emission by par- 
ticles with different energies may produce a continuous spec- 
trum, where the spectral lines are barely visible or completely 
negligible. 

The approximation of the cyclo-synchrotron emission pro- 
vided by Marcowith & Malzac ( 2003 1 has been compared with 
the results of the precise numerical computations showing dis- 
crepancies that are less than 20%. Of course, the discrepancy 
between our continuous approximation and the precise calcu- 
lation of (discrete) harmonics can be very large, if we com- 
pare our formula with the emission level between two well- 
separated harmonics. On the other hand, our main goal is to 
derive a formula that always provides the correct value of the 
total emitted energy. This is achieved through the normalization 
term (Eq. that independent of the values of the constants a 
& b, always provides the correct cooling ratio. This construc- 
tion of the formula introduces a freedom in manipulating of 
the spectral shape. Therefore, by choosing appropriate values 
for the parameters a and b we can approximate the spectra at 
different energies rather well. 

In Fig.^we compare our new formula with the old relation 
and the approximation of the first ten harmonics provided by 
Marcowith & Malzac (■2003J . A few general conclusions can 
be drawn from this comparison: 
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Fig. 1. The comparison between different approximations of 
the cyclo-synchrotron power spectrum of a single particle at 
four different energies. The thin solid line in each panel shows 
the polynomial approximation of the first ten harmonics pro- 
vided by Marcowith & Malzac ( I2003I I. The continuous approx- 
imation derived in this work is shown by the bold solid line. 



The simple approximation proposed by Ghisellini et al. ( I1998> 
is shown by the thin dashed line. Note that the last approxi- 
mation gives the correct cooling rate only if integrated from 
frequency vb (indicated by the vertical dotted fine). 
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• For very low particle energies (y6 <K 0. 1), most of the power 
is emitted in the first harmonic, i.e. in a very narrow fre- 
quency range. For this energy range, our approximation 
shows the biggest disagreement in comparison to the ap- 
proximation of the first harmonic. However, in compared 
to the old formula, our approach gives significantly better 
results. 

• In the intermediate energy range (0.1 < /3 < 0.5) the con- 
tribution of high-order harmonics to the total emission be- 
comes more important. Moreover, with increasing particle 
energy, the emission from each harmonic spreads over a 
wider frequency range. The spectrum transforms from a set 
of discrete lines into a continuous emission. Therefore, our 
approximation as a continuous function becomes more and 
more accurate with increasing particle energy. 

• Finally, in the range of relatively high energy (0.5 < /3 < 
0.9), the emission is dominated by the high-order harmon- 
ics. Our approximation gives the best agreement with the 
approximation for the first ten harmonics. The old formula 
also provides a relatively good approximation in this en- 
ergy range. However, the total radiated power calculated 
from this formula is correct only if integrated starting from 
frequency vb- Therefore, it cannot describe the spectrum 
below Vb, since it overestimates the total emitted energy if 
integrated from v^in The approximation of the first ten har- 
monics does not provide a correct value of the total emitted 
energy, either due to the fact that the ten harmonics consid- 
ered provide a good description only up to v ~ 4-vb, while 
there is still considerable power above, too. For this particu- 
lar particle energy range, many more harmonics are needed 
to provide the correct value of the total emitted energy. 

Above jS - 0.9, the synchrotron power spectrum can be used; 
however, this coefficient also requires some correction, which 
we describe in the next section. 

3. Correction of the synchrotron power spectrum 

The synchrotron power spectrum from a single particle in a 
random magnetic field, integrated over an isotropic distribution 
of pitch angles, has been derived by Crusius & Schlickeiser 
SWH6» and GGS88: 
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where x - v/(3y^VB) and Ky{x) is the modified Bessel function 
of order y. This formula does not provide the correct cooling 
rate for y < 15, where it overestimates the total emitted en- 
ergy. Since the spectral shape, however, is approximately cor- 
rect even for low particle energies (y < 15), we multiply the 
original formula by a simple correction term that only depends 
on particle energy 
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Fig. 2. The approximated cyclo-synchrotron power spectrum 
and the corrected synchrotron spectrum for p = 0.9. The thin 
solid line shows the approximation of the first ten harmonics 
provided by Marcowith & Malzac (2003 1, the continuous ap- 
proximation derived in this paper is shown by the bold line, the 
synchrotron power spectrum derived by Crusius & Schlickeiser 
(>1986 ) and GGS88 is shown by the short dashed line and the 
long dashed line shows the corrected synchrotron spectrum. 

The difference between the cooling rate iyc) and the frequency 
integrated synchrotron power disappears for y » 15; therefore, 
the correction term reduces to unity for high particle energies. 

In Fig. 121 we compare our approximated cyclo-synchrotron 
power spectrum with the corrected synchrotron spectrum at 
P - 0.9. Since the correction term for the synchrotron formula 
only depends on the particle energy, the correction only affects 
the normalization of the synchrotron spectrum. The figure also 
shows, that our approximation joins relatively smoothly to the 
corrected synchrotron spectrum at v » v^. However, in order to 
achieve this smooth connection, we had to modify the constant 
a (Eqs|2and0, that controls the spectrum shape at v » vb- 
The constant h in our formula controls the spectral shape for 
P < 0.5. Note that any modification of the parameters a or b 
changes the spectral shape and thus, in principle, also the total 
emitted energy. This problem has been solved through the nor- 
malization term (Eq. 0, which also contains the parameters a 
& b, and thus changes the level of the spectrum in order to keep 
the correct value of the total emission. 

4. Particle heating through the absorption of 
photons 

We present a simple application for our approximation of 
the cyclo-synchrotron power spectrum and the corrected syn- 
chrotron emission coefficient. We analyze the amount of en- 
ergy gain corresponding to the cyclo-synchrotron absorption 
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process. This process may lead to an efficient exchange of 
the energy, and may therefore provide a very powerful heat- 
ing mechanism for low energy particles. This kind of heating is 
a stochastic process. This process of competing radiative cool- 
ing and radiative heating, through emission and re-absorption, 
respectively, leads to the accumulation of most particles around 
the equilibrium energy (ye), where heating and cooling are in 
balance. In an ideal case, such a competition would transform 
any initial particle distribution into a thermal Maxwellian dis- 
tribution with its maximum at the equilibrium energy. 

In our test, we investigated the self-absorption of the ra- 
diation field produced by electrons with a power law energy 
distribution 



N(y) oc y-". 



(11) 



located inside a homogeneous spherical volume that is filled by 
a tangled magnetic field. We compare the efficiency of the heat- 
ing and the cooling processes and discuss the possible values 
of the equilibrium energy for different slopes of the electron 
energy distribution. 

The emission coefficient for any electron energy distribu- 
tion N(y) is defined by 



j(y) 



4;r/ 



N{y)P{v,y)dy, 



(12) 



where P(v, y) is either the cyclo-synchrotron or the corrected 
synchrotron power, depending on the value of y. The self- 
absorption coefficient for any particle distribution is given by 
LeRoux(1961i 



k(v) 



N(y) d 
yp dy 



[ypP(v, y)] dy. 



(13) 



For the calculations we assume a constant intensity of the 
radiation field inside the source, which is approximated by 



/(v) = C^{l-exp[-^(y)7?]}, 
k(v) 



(14) 



where R is the source radius. The value of the parameter C de- 
pends in principle on the position within the source. In the cen- 
ter of our homogeneous source C = 1 should be used, whereas 
on the source surface we should apply C - 1/2 (Gould . 1979 1. 
However, for the sake of simplicity we assume a constant value 
of the intensity for the whole absorbing region, and use an av- 
erage value of C = 3/4. 

Finally, the absorption or heating efficiency for any radia- 
tion field is described by (GhiselUni & Svensson ll991> as 



1 1 d 

nieC^ yp dy 
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(15) 



We can easily estimate the absorption efficiency in an asymp- 
totic case where the particle energy is characterized by y — > 1 
(/3 <K 0.1). In this case, most of the energy is emitted in the first 
harmonic, i.e. in a very narrow frequency range (see Fig.^for 
/3 - 0.01). Therefore, the cyclo-synchrotron power spectrum 
can be approximated by a 5-function 



which results in a constant value of the absorption efficiency 

4 /(v^) / p2\ 

iair ^ 1) = :;co-tUb- j\3y + — | ^ constant. (17) 
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Pc(y, r -> 1) = meC ydp) 6{v - vb) 
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Comparing this result with the cooling efficiency, which is 
always proportional to p^ (Eq. |3}, we see that /or low energy 
particles, heating due to the self-absorption will always over- 
come radiative losses. In other words, for a power law particle 
energy distribution, the heating term is proportional to p^ for 
high values of p (in the self-absorbed regime); on the other 
hand, in the low energy limit /:> ^ (or y — > 1), the heat- 
ing term must always be constant. This implies, that there al- 
ways is only one specific energy value for which the heating 
and cooling terms are equal. To illustrate this point, let us com- 
pare this result with the equilibrium solution {N{y) oc -y^^] pro- 
posed in the "turbulent reactor" scenario. In this case, the (an- 
alytical) solution is obtained by considering only the relativis- 
tic regime and further assuming an infinite source (i.e. infinite 
self-absorption frequency). In other words, the radiation inten- 
sity is assumed to be proportional to v'l^ for all frequencies, 
with a normalization that depends on the slope of the electron 
distribution (see also Rees 1967 and Mc Cray 1969 1. If one as- 
sumes that the particle distribution is a power law, but truncated 
at some low energy y^i^ (to self consistently use the emissivity 
and absorption processes in the relativistic regime) then the ra- 
diation field is not oc y^l^ in the entire frequency range, but only 
above ~ y^ . v^. Below this frequency, I(v) oc (e.g. Rybicki 

6 Lightman 1979 1, making the electrons gain energy (through 
absorption) at a slightly higher rate than what is found by as- 
suming /(v) oc v^^'^. If instead one assumes a particle distribu- 
tion extending towards mildly relativistic energies, then one ob- 
tains that 7fl, for small particle energies becomes constant and 
thus clearly larger than the extrapolation of the % oc law. 
Concerning the other energy extreme, a source of finite size 
becomes transparent at some finite value of the self-absorption 
frequency v,. Therefore the radiation field inside the source for 
frequencies close to and above v, is no longer proportional to 
v^^^, with a corresponding decrease in the energy gain rate for 
high-energy particles. These are the reasons a particle distribu- 
tion A^('>') oc y^^ is not an equilibrium solution. For any given 
value of the slope of the particle distribution, equilibrium is al- 
ways achieved at only one specific energy. 

In Fig. |3l we compare the absorption efficiency with the 
cooling rate of a power law particle energy distribution from 
ySmin = 0.01 to y8„,x = 0.9999 (r^ax ^ 70) with three different 
slopes (n = 2, n = 3, and n - 4). We performed the computa- 
tions using three different approaches for the calculation of the 
cyclo-synchrotron power spectrum in order to test the formulae 
derived in this work. 

First, we used the old formula for the cyclo-synchrotron 
power spectrum (Eq. |2} and the standard, uncorrected syn- 
chrotron emission coefficient (Eq.|S|l. The result of these com- 
putations is presented in the upper left panel of Fig.|3l In this 
particular case, the transition from the old cyclo-synchrotron 
power spectrum to the uncorrected synchrotron emission at 

7 = 2, produces a clearly visible discontinuity in the absorption 
efficiency. The discontinuity is related to the fact, that the un- 
corrected synchrotron emission overestimates the total emitted 
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Fig. 3. The efficiency of absorption (solid lines) versus the cooling efficiency (long dashed lines) for the power law particle energy 
distribution ifiram = 0.1 — » ySmax = 0.9999, j ^ 70) for n - 2, n - \, and n = 3 (the last only in the main panel - c). The left 
upper panel (a) shows the absorption efficiency calculated from the approximation of the cyclo-synchrotron emission derived 
by Ghisellini et al. f 1998^ and the standard synchrotron emission coefficient provided by Crusius & Schlickeiser (I1986> and 
GGS98. In the upper right panel (b), we show the absorption efficiency calculated only from the corrected synchrotron emission. 
The dash-dot lines on this panel indicate the result obtained from uncorrected synchrotron emission. In the main panel (c). the 
new cyclo-synchrotron approximation and the corrected synchrotron emission were used. The short dashed lines in each panel 
indicate a constant value of the absorption efficiency estimated for y ^ 1 (Eg.llTL The equilibrium energies are indicated by 
dots. Note that in the main panel we also show the absorption efficiency calculated for ySmi,, = 0.9 and n = 3 (bold dotted Une). 



energy for y < 15. Moreover, for low energy particles (/? < 0. 1) 
the level of the absorption efficiency does not agree with the 
constant level estimated from the ^-approximation (Eq. I17> . 
which is indicated by the horizontal lines. 

In the second test, presented on the main panel of Fig.|3 we 
used our new approximation of the cyclo-synchrotron power 



spectrum and the corrected synchrotron emission coefficient. 
Since both expressions provide the correct value for the total 
emitted energy, the transition from one spectrum to the other 
(this time at /? = 0.9) does not produce any discontinuity in the 
absorption efficiency. Moreover, with our new approximation 
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of the cyclo-synchrotron emission, the absorption efficiency at 
low particle energies agrees well with the (5-approximation. 

Finally, we used only the corrected synchrotron emission 
coefficient for the whole energy range (right panel in Fig. 
No discontinuity is present, since we used only one formula. 
However, this approach, as well as the old cyclo-synchrotron 
power spectrum in the first test, do not agree very well with the 
(^-approximation. Note that for low energy particles, the uncor- 
rected synchrotron emission gives an absorption efficiency that 
is a few orders of magnitude larger than the efficiency obtained 
from the other approximations. 

These three approaches for the calculating of the absorption 
efficiency, qualitatively give the same results, but differ in the 
quantitative details, which indicate that the new approximation 
for the cyclo-synchrotron power spectrum together with the 
corrected synchrotron emission coefficient, provides the most 
precise description. 

Our tests indeed show that the absorption efficiency for the 
very low-energy particles becomes independent of the particle 
energy and is significantly higher than the cyclo-synchrotron 
cooling ratio. Therefore, self-absorption will always cause a 
strong heating of the low-energy particles. The equilibrium en- 
ergy depends strongly on the slope of the particle spectrum. 
For n - 2 the equilibrium energy is very close to the peak in 
the absorption efficiency (the jaip) curve), which is related to 
the maximum in the self-absorbed spectrum (v,). If the particle 
spectrum is steeper (n - 4), the equilibrium is taken at a lower 
energy, but close to the energy of those particles emitting at the 
peak of the synchrotron spectrum. The n = 3 case is particu- 
larly interesting, since it corresponds to the previously claimed 
equilibrium solution. Contrary to this claim, heating and cool- 
ing also balance in this case only at a specific energy. This is, 
on one hand, due to the finite size of the source that limits the 
range of possible momenta of electrons that emit and absorb 
radiation efficiently and on the other hand, more importantly, 
due to the trans- and sub-relativistic regime where low energy 
particles always gain more energy than they loose. 

The value of the equilibrium energy strongly depends on 
the minimum and maximum energy of the particles. In our 
tests the equilibrium energy for n < 3 depends on the self- 
absorption frequency (v,). However, for relatively low value of 
Tmax (e-g- 7 instead of 70 in our particular calculations), the 
emission should be absorbed at all frequencies of a completely 
optically thick source. In such a case, the equilibrium energy 
depends directly on ymax- For n > 3 the equilibrium energy de- 
pends on the minimal energy of the particles. In Fig.|3we show 
the heating efficiency calculated for a relatively high value of 
the minimum particle energy ySmin = 0.9. The value of the equi- 
librium energy increases with the increasing minimum energy. 

Note that for sake of simplicity, our tests assume a sta- 
tionary state, where the equilibrium energy is simply given by 
the equilibrium between the heating and cooling rates. In real- 
ity, the system evolves and the physical conditions inside the 
source change. The initial power law, or any other particle dis- 
tribution, will be transformed into a thermal or quasi-thermal 
spectrum (see GGS88). Also, the equilibrium energy in such 
an evolving source may be different from the energy estimated 
from the simple stationary analysis. This does not change that 
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there is only one preferred equilibrium energy around which 
most of the particles will be accumulated, forming thermal or 
quasi-thermal distribution. A complete description of this time 
dependent evolution will be the main focus of our future study. 

5. The self-absorbed spectrum 

The detailed analysis of the cyclo-synchrotron emission shows 
that the self-absorbed part of the spectrum generated by the 
electrons with a power-law energy distribution, is slightly dif- 
ferent from the well known power law relation /(v) oc v^l^. 
In this section we analyze the emission of our homogeneous 
source, assuming different slopes of the particle spectrum, and 
discuss the reasons for the deviations from the standard v^^^ 
spectrum. 

The observed intensity of the emission from the homoge- 
neous spherical source is given by 

where t = 2Rk is the optical depth (e.g. Bloom & Marscher 
1996 1. We calculate the intensity for a range of power law 
slopes, starting from « = -1 up to n = 25/3 with finite steps of 
An = 2/3. Our results are shown in Fig.|4] For n < -2/3, the 
spectral index is constant at or = -2, in almost the whole self- 
absorbed part of the spectrum. In the range -2/3 < n < 5/3, 
the slope of the spectrum changes from a = -2 to a - -I for 
V < vb, and from a = -2 to o- = -2.5 for v > vb- Above the 
limiting value n = 5/3, the spectral index remains constant 

a — —I for V < Vb, 
a - -5/2 for v > vb- 

Note that, for relatively steep particle distributions (n > 17/3 
in this particular case) close to the low frequency cut-off, the 
spectrum may go over into a power law with the index a - -n. 
Therefore, in such a case the self-absorbed part of the spectrum 
is described by three different indices a = -n — » -1 — » -5/2. 
It is relatively easy to understand this specific evolution of the 
self-absorbed spectrum, if we analyze the limiting cases n < 
-2/3 andn > 5/3. 

In the first case (« < -2/3), the emission in the whole fre- 
quency range is dominated by the synchrotron radiation pro- 
duced by the highest energy electrons. The synchrotron power 
spectrum of a single particle or monoenergetic population of 
the particles can be approximated below the peak frequency 
(vp(y)) by 

P,iv,y)ocv^l^ y-^l^. (19) 

Integrating this formula over the power-law electron energy 
spectrum with the index n < -2/3, we obtain 

j{v) oc v'^^ for Vmin(rmax) « V <K Vp(7'niax)- (20) 

Calculating the absorption in the same way, Eq.^] we obtain 

k{v) oc v'^l^ for V,ni„(yn,ax) <^ V <^ Vpiy^^^), (21) 

and this gives /(v) = j{v)/k{v) oc in the self-absorbed part 
of the spectrum. 
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Fig. 4. The self-absorbed part of the cyclo-synchrotron and the 
synchrotron emission from a homogeneous spherical source. 
The intensities were calculated for a range of different indices 
of the power law particle energy spectrum, starting from y' up 
to y"^^/^ with a finite step size An = 2/3. The spectrum on 
the top [/(v) oc in almost the entire frequency range shown] 
corresponds to n = -1, whereas the lowest spectrum was cal- 
culated for n -25/3. The limiting values of the spectral indices 
presented in this figure are discussed in the Sect.|5] 

In the second limiting case (« > 5/3), the spectral index 
above vb is equivalent to the well-known solution (a = -5/2), 
and we only discuss the reason for the flattening of the spec- 
trum (a - -1) below this frequency. Around vb the emis- 
sion is dominated by the cyclo-synchrotron radiation of the 
low-energy particles. However, for v vb the emission be- 
comes dominated by the tail of the synchrotron emission of 
the high-energy particles. Therefore, we can again use the ap- 
proximation of the single-particle power spectrum (Eq. I19> . 
Integrating this approximation over the power law electron 
spectrum and neglecting the lower integration boundary, we ob- 
tain j(v) oc v'^-^y^(,x According to Eq.|5lthe maximum energy 
is directly related to a given frequency v ~ l/ymax- Therefore, 
we obtain 

j(v) oc v" for Vmin(7tnax) « V « Vfi. (22) 

In the same way, we can integrate the absorption coefficient 
Ea.ll3lobtaining k ~ v"^^"'')'max^^''^ which gives 

^(V) oc v"-l for Vmi„(rmax) « V <SC Vg. (23) 

Finally, this gives /(v) = j(v)/k(v) oc in the self-absorbed 
part of the spectrum below the frequency vb- 

The low-frequency emission (v <K vb) can be absorbed 
only by the high-energy particles. For relatively steep parti- 
cle spectra, the density of these particles can be too small to 
efficiently absorb the low-frequency radiation. Therefore, the 



source may again become optically thin and the spectral index 
equivalent to the index of the emission coefficient /(v) oc v" 
(Eq.EU. 

When calculating the spectra presented in Fig.|3 we used 
our new formula for the cyclo-synchrotron power spectrum 
and the corrected synchrotron-emission coefficient. There is 
no smooth transition between these formulae at /3 = 0.9 for 
y Vb (see Fig. but this is barely visible in our spectra. 
The effects discussed in this section depend mostly on the syn- 
chrotron emission; therefore, the corrected synchrotron emis- 
sion can be used in whole energy range to get almost identical 
results. 

The modifications of the self-absorbed spectrum that we 
discuss appear at relatively low frequencies. Therefore, in most 
astrophysical objects such effects are not observable. However, 
some effects might be visible in some specific physical condi- 
tions. One example might be the synchrotron radiation (in its 
self-absorbed portion) produced by steep power-law distribu- 
tions of particles in highly magnetized sources {B ~ 10^^'" 
[G]). Note that isotropic distribution of the pitch-angles, as- 
sumed in order to derive the emissivity formulae presented in 
this paper, may not always be valid, especially in very highly 
magnetized sources. 

6. Discussion 

We have derived a new approximation for the cyclo- 
synchrotron power spectrum of a single particle and com- 
pared it with the approximation to ten first harmonics of the 
cyclo-synchrotron emission provided recently by Marcowith 
& Malzac (2003). In comparison to the other approaches, our 
approximation self-consistently provides the correct value of 
the total emitted energy over the whole range of the particle 
energies. Moreover, our approach describes the spectrum of the 
emission in the range 0.1 < /3 < 0.9 relatively well. Finally, the 
approximation provides a relatively smooth connection with 
the corrected synchrotron emission at v > for /3 - 0.9. 

All these results are useful when one needs fast computa- 
tional tools to derive the cyclo-synchrotron emission and ab- 
sorption, instead of using the exact expressions, which require 
much more computing time. 

The application we will pursue is to study in detail the 
evolution of the emitting particle distribution subject to ac- 
celeration and/or injection of new particles, radiative and 
Coulomb cooling, and heating due to the synchrotron self- 
absorption. For the moment, we have instead analyzed a sim- 
pler process that, however, requires a careful treatment of the 
trans-relativistic regime. We have demonstrated that a power- 
law distribution of electrons emitting and absorbing cyclo- 
synchrotron photons can never be a steady solution. Energy 
losses and gains always equal each other at a particular energy 
-yc, and not over a range of energies. This is contrary to previ- 
ous claims that a distribution N(y) oc can be an equilibrium 
solution (e.g. Kaplan & Tsytovich 1973). The reason is that, 
at low enough, sub- or trans-relativistic energies, energy gains 
always exceed losses and, at the other extreme, the absorption 
becomes less efficient because the source becomes transparent 
(unless it is infinite in size). Particles therefore will tend to ac- 
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cumulate at jg, changing the shape of the energy distribution 
they initially belonged to. 

For the steep particle spectrum (n > 3), the equilibrium 
energy (y^.) is small. This motivated our detailed study of the 
trans-relativistic cyclo-synchrotron regime. However, we have 
shown that different approximations lead to very similar re- 
sults for the amount of energy gains experienced by the par- 
ticle. This is due to the fact that in calculating this quantity 
we must consider frequency-integrated expressions, with the 
consequent loss of details concerning the shape of the power 
spectrum. What matters is mainly that the frequency integrated 
spectrum equals the correct cooling rate (i.e. oc p^). 

Exchanging photons through emission and absorption al- 
lows particles to exchange energy, independently of Coulomb 
collisions. This is a very important thermalization process in 
those magnetized, hot and rarefied plasma where Coulomb col- 
lisions are rare. Cyclo-synchrotron absorption transforms an 
initially non-thermal distribution into a Maxwellian in just a 
few cooUng times (GGS88). 
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